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Abstract: We use a variational technique to study heavy glueballs on gauge configura- 
tions generated with 2+1 flavours of ASQTAD improved staggered fermions. The varia- 
tional technique includes glueball scattering states. The measurements were made using 
2150 configurations at 0.092 fm with a pion mass of 360 MeV. We report masses for 10 
glueball states. We discuss the prospects for unquenched lattice QCD calculations of the 
oddballs. 
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1 Introduction and motivation 

Nothing is more symbolic of the difficulty of solving QCD, than the fact that, while glueballs 
are central to the understanding of non-perturbative QCD, there is currently no definite 
experimental evidence for their existence. After much work the glueball spectrum [1] 
in quenched QCD was mapped out by Morningstar and Peardon and collaborators [2, 3]. 
Results for higher spin glueballs have been reported by Meyer and Teper [4, 5] (see also [6]). 

There has been much less work done on studying the effect of sea quarks on the glueball 
masses [7-10]. Glueball calculations suffer from a severe problem with the signal to noise 
ratio that requires high statistics. Also, some of the more sophisticated algorithms used 
in quenched QCD calculations that improve the signal to noise error, do not work for 
unquenched calculations [5, 11]. 

The properties of glueballs can be elucidated by studying the experimental decay and 
production of flavour singlet mesons. Some of these analyses find that the glueball degrees of 
freedom have similar masses to the quenched glueballs [12], while others find very different 
masses [13, 14]. For example, one analysis [13] of the decay properties of the h states 
suggested that large unquenching effects moved the quenched h glueball from 2.6 GeV 
to 1.4(1) GeV, close to the experimental mass of the 7/(1405) meson. 

Some have argued (e.g. [15]) that there is an additional state over the two expected 
qq mesons in the ++ flavour singlet mesons: /o(1370), /o(1500), and /o(1710). The mass 
of the ++ glueball inferred from quenched QCD [2] is 1730(50) (80) MeV. This could be 
important in understanding the above mesons if one assumes that unquenching effects 
m the ++ glueball are small. Other groups, e.g. [14], have argued that unquenching of 
the ++ may be large and that the mass in QCD is close to the a. There is also some 
controversy as to whether the the /o(1370) is a real meson state [16], although it is listed in 
the Particle Data Group (PDG) summary tables. Lattice calculations will eventually have 
to deal with the decays of the /o mesons and possible coupling of the /o(980) and /o(600). 
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The mixing of the glueball degrees of freedom with flavour singlet mesons has meant 
that many lattice groups [17-21] have been studying the p meson and P-wave states such 
as the ai(1260), b\ (1235) meson using a variety of techniques so as to understand how to 
deal with resonances in lattice QCD calculations. 

So, in summary, there are no hadrons where glueball degrees of freedom have been con- 
firmed. Looking to the future, there are ongoing experiments that are searching for glueball 
degrees of freedom, such as the BES III experiment [22]. The BES III experiment has pre- 
liminary [23] results for two new states with masses close to the 0^ and 1 H quenched 
glueball masses. In 2018, the PANDA experiment [24] will search for heavy glueballs with 
masses under 5.4 GeV. In particular they will look for oddball glueballs with exotic J PC 
quantum numbers (0 -1 , 2 H , 1 , and 3 h ) which are not allowed in quark models 
for quark-antiquark mesons. In a quenched QCD calculation, Morningstar and Peardon [2] 
found two glueballs with exotic J PC = 2 + ~ and + ~ with masses 4140(50) (200) MeV and 
4740(70) (230) MeV respectively. As we review later, other quenched glueball studies have 
not seen these states [25]. 

One motivation for studying glueballs heavier than 3 GeV is that there could be reduced 
mixing with flavour singlet quark states. There have been speculations that there are no 
further, or a reduced number of, light meson states above 3.1 GeV [26, 27], based on string 
breaking of the heavy quark potential. Swanson [28] critiques the use of screened potentials 
in hadron spectroscopy. The heaviest light meson in the PDG summary table [29] is the 
/e(2510) with a mass of 2.469 GeV. It is possible that there may be no experiments capable 
of producing light mesons beyond 2.5 GeV. It could be that heavier mesons made from light 
quarks have large widths, so it becomes difficult to extract the masses from experiment. 
Even though the signal to noise ratio is worse for heavier glueballs than for light glueballs, it 
would simplify lattice QCD calculations if heavier glueballs do not mix with quark degrees 
of freedom, because it would be simpler to identify a pure glueball state. For glueball masses 
above 3 GeV, there is the possibility that the glueballs may mix with charmonium states. 
However, Page [27] suggests that the mixing between charmonium states and glueballs will 
be small, because the charm quark is heavy and hence it is difficult to excite a charm loop. 

The hadron spectrum collaboration [30] has computed the spectrum of excited mesons 
and baryons, although typically at a single lattice spacing and heavy quark masses. They 
don't report any meson masses above 2.8 GeV. Studying such masses is already an impres- 
sive achievement, and it is not clear whether they can investigate the meson spectrum at 
even higher masses. 

Hagedorn conjectured that the density of light hadrons goes like e m l T where m is the 
mass of the light hadrons and T is a constant [31, 32]. Cohen and Krejcirik have recently 
reviewed [33] the evidence that the number of light hadrons agrees with the Hagedorn 
conjecture. 

There may be indirect ways of determining whether there exist light mesons in the 
regime 3 to 6 Gev, using information from thermodynamic studies. For example, the hadron 
resonance gas model (HRGM) is used in the phenomenology of heavy ion experiments 
and has successfully reproduced some results from lattice QCD calculations at non-zero 
temperature. The HRGM depends on the number of mesons and baryon states. Typically 
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p 


L 3 x T 


ami 


ami 


N cfg 


Ntraj 


r /a 


7.095 


32 3 x 64 


0.00775 


0.031 


2150 


12900 


5.059(10) 



Table 1. Summary of the parameters of the ensemble used in this calculation. We use the same 
convention for the quark masses as used by the MILC collaboration. The configurations used are 
taken every 6 trajectories. 

the meson and baryon states listed in the PDG are used, however Chatterjee et al. [34] 
have studied a different density of states, with mesons with masses higher than those listed 
in the PDG. Majumder and Muller use the results from lattice QCD [35, 36] calculations 
with the HRGM to claim the existence of new light resonances [37]. Megias et al. [38] 
study the effects of the additional hadrons predicted by the quark model on the results of 
the HRGM. 

The PANDA experiment is performing some Monte Carlo studies of glueball production 
using glueball decay widths of 10 MeV [39]. It will be important to check some of these 
speculations about the density of light mesons and the mixing of glueball and charmed 
hadrons. 

Oddball glueballs will not mix with light mesons with non-exotic J PC quantum num- 
bers. Even if the heavier glueball states do not mix with quark degrees of freedom, there 
is the possibility that the glueballs decay to other glueballs. The calculation reported here 
explicitly includes glueball scattering states to check for this. 

In [10], the results for ++ , and 2 ++ were presented on the gauge configurations 
also used in this study 1 . A powerful variational code was developed in [40] that included 
two glueball scattering states and was used to study glueballs in SU(N) quenched QCD. 
In this paper, we use this code to extend the study in [10] to look at glueballs with new 
J . We will mostly study the glueballs with J PC quantum numbers not studied in the 
earlier work [10]. 

There are a number of reviews on the theory of glueballs [41, 42], lattice QCD calcu- 
lations of glueballs [43], and experimental searches [44] for glueball degrees of freedom. 

2 Details of the lattice QCD calculation 

The unquenched lattice QCD calculations used the ASQTAD improved staggered fermion 
action [45, 46] and the one link Symanzik improved gauge action [47]. These are the same 
actions used by the MILC collaboration [48]. The configurations were described in one 
paper on glueball masses [10] and one on the masses of the rj and rf mesons [49]. The 
parameters of the calculation are shown in Tab. 1. The lattice spacing a was determined 
from the r§ja parameter calculated from the heavy quark potential and the value of ro 
(0.4661 fm) determined by the HPQCD collaboration [50] on the configurations generated 
by the MILC collaboration [48]. 

Glueball masses were extracted with a variational calculation. The variational tech- 
nique used in this work was described in [40]. Here we provide a brief summary. The 

In particular, this work uses the set of configurations referred to as the fine ensemble in [10]. 
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basis states for the variational procedure are made from operators that transform under an 
irreducible representation of the cubic lattice group and have specified parity and charge 
conjugation quantum numbers. The irreducible representations of the cubic group, which 
is the discrete rotational symmetry group of the lattice, are conventionally called Ai, A2, 
E, T\ and T2. Representations A\ and A2 have dimension 1, E has dimension 2 and T\ and 
T2 have dimension 3. Below, we show how the continuum spin can be obtained from mass 
eigenstates transforming according to an irreducible representation of the cubic group. 

On a given timeslice, a single glueball operator with well-defined rotational quan- 
tum number is a prescribed linear combination of traced Wilson loops of a given shape. 
Eigenstates of parity are obtained by considering the starting linear combination and the 
reflected one, while eigenstates of charge conjugations are given by the real part (C = 1) 
and the imaginary part (C = —1) of that combination. In our calculations, we considered 
operators containing shapes from length 4 to length 10. In addition to single glueball oper- 
ators, the basis states include two-glueball scattering states (which involve products of two 
basic shapes) and bi-torelon operators (products of two loops winding in a compact spatial 
direction; note that these operators are not necessarily straight lines). In order to keep the 
number of operators under control, we limited the number of single glueball operators to at 
most 30. The shapes we included in our calculations are those that contribute in as many 
channels as possible and, barring shapes that for symmetry reasons do not contribute to 
given channels, the same basis shapes have been used for all the quantum number assign- 
ments. Other constraints imposed on the basis operators were that the variational basis 
must contain shapes of all lengths and that a given shape can contribute at most once in a 
given channel. This allowed us to include a broad range of shapes and lengths in our basis. 
A Mathematica program systematically constructed the basis operators with the requested 
constraints. A summary of the number of operators we have used for each quantum num- 
ber assignment is reported in Fig. 1. Additional operators were generated by an improved 
blocking algorithm [51]. We used 6 different blocking and smearing levels, which implies 
that for each channel our basis is 6 times larger than the number of operators quoted in 
Fig. 1. 

Consider two basis operator $ a and In the lattice QCD calculations the correlators 
between operators separated by time t are computed and define the correlation matrix 

Cafi(t) = 52(0\& a (t + T)* P (T)\0) • (2.1) 

T 

The masses and amplitudes of glueball states can be extracted from the generalized 
eigenvalue problem on the measured correlation matrix: the optimal operators (i.e. those 
that almost coincide with pure states) are the eingenvectors of 

C7(t) = C-\0)C{i) . (2.2) 

These interpolating operators are a linear sum of the basis vectors: 

**(*) = • (2-3) 

a 
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Figure 1. Number of basis operators (before blocking and smearing) in each channel, for single 
glueballs, scattering states and bi-torelon states. 



The mass of the state is extracted by cosh fits of the correlation matrix in the optimal 
basis: 

C u (t) = \a\ 2 cosh {nut - N T /2) , (2.4) 

where Nt is the length of the lattice in the time direction and the cosh functional form is 
a consequence of the usual exponential decay in a lattice with periodic boundary condition 
in the time direction. 

In general, glueball correlators are very noisy and this limits the usefulness of numerical 
correlators to short time separations. However, although Eq. (2.4) is only valid at large 
t, if the overlap with an Hamiltonian state is almost perfect, it is possible to extract a 
reliable value for the mass at short time separation, since the decay is largely dominated 
by a single state. For this to be true, a careful construction of the variational basis is 
paramount. Whether an optimal state is a good approximation of the Hamiltonian 
eigenstates can be checked by looking at the value of the overlap |cj| 2 : the closer this number 
is to one from below (with one being the unitarity limit), the better is the variational 
calculation. In addition, it is important to estimate the contributions to the mass coming 
from scattering and torelon states. For the eigenstate corresponding to a true spectral mass, 
these contaminations must be absent. The contribution of scattering and torelon states 
for a given optimal state can be resolved by looking at the relative length of its projection 
onto the space spanned respectively by the scattering and bi-torelon basis operators. A 
summary of the resulting relative projections is shown in Fig. 2. For the lightest states 
extracted in each symmetry channel we show the relative projection onto the different types 
of basis operators (single glueballs, scattering and bi-torelon operators). 

We do not include any fermionic scattering states in this calculation. Fu has recently 
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Figure 2. Summary of the projections relative to the subspace spanned by different types of basis 
operators for the lightest states in the spectrum. 



reported [52] preliminary results for the decay width of the ++ a into tttt. New techniques 
for calculating the masses of scalar mesons using lattice QCD have been developed [53, 54]. 

The results come from bin sizes of 50 configurations. We have checked for autocorrela- 
tion effects by choosing different binning blocks down to 5 configs. We found no evidence 
for correlations. 



3 Results 

We reduced the variational correlators Eq. (2.1) truncating the diagonal matrix Eq. (2.2) to 
the 5 highest eigenvalues (corresponding to the 5 lowest masses). The diagonal correlators 
are fitted to the fit model in Eq. (2.4). The functional form is a good approximation if the 
overlap |q| 2 ~ 0(1) and the x 2 /dof ~ 0(1)- All the states we were able to measure have 
0.65 < |cj| 2 < 1.1, with order 0.1 errors and < x 2 /dof < 1.25. 

The correlator for the lightest glueball (the channel, which in the continuum 

encodes the scalar ++ glueball) and its first excitation is shown in Fig. 3 on a logarithmic 
scale. The effective mass plateaux of these two states with the corresponding estimate of 
the fitted masses from the fit in Eq. (2.4) are reported in Fig. 4. A similar plot for the E ++ 
channel corresponding to the continuum tensor glueball is shown in Fig. 5. For the At , we 
have investigated alternative variational calculations that include only some of the subsets 
of our operators (only single glueball operators or only scattering and torellon operators). 
We note that in both cases the mass of the ground state is the same within errors. This 
phenomenon is different from what happens in pure gauge, where for sets excluding single 
glueball operators the ground state mass is roughly twice the mass of the ground state one 
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Figure 3. Normalized correlator in log scale for the two lightest states in the spectrum of the +H 
channel. 
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Figure 4. Effective mass plateaux for ground state (circles) and the first excitation (squares) in the 
++ channel. The fitted values are shown with solid lines (the dashed lines represent one standard 
deviation) . 



finds when single glueball are included in the variational basis [40]. This might be due to 
the explicit breaking of centre symmetry in the presence of dynamical fermions. Although 
the coupling of glueball states and single torellon states was found to be negligible in [8], 
it is possible that it becomes more relevant closer to the chiral limit. To investigate this 
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Figure 5. Normalized correlator in log scale for the lightest state in the spectrum of the E ++ 
channel. 

coupling systematically would require variational calculations with a basis including single 
torellon operators at various quark masses close to the chiral limit. A similar calculation 
is beyond the scope of this paper. 

In Fig. 6 we plot the glueball masses for each lattice representation. A summary of 
the spectrum is reported in Tab. 2. We define glueball states as those with an overlap of 
at least 65% onto the single glueball operators. 

We also include the results from Chen et al. [3] from a quenched lattice QCD calculation in 
the continuum limit. Fig. 6 does not show any significant unquenching effects, although the 
continuum limit of the unquenched results should be taken before a definitive statement 
can be made. 

4 Classifying the glueballs into continuum angular momentum states 

Each of the states comes from an operator that is a linear combination of ones in the 
variational basis. The projection on the 3 subsets can be used to identify single glueball 
states from scattering or bi-torelon states. Only the , and channels had 

O(40)% overlap with the scattering subset of operators. Also only the A~[ + and 
channels had contaminations of O(40) % overlap with bi-torelons. We find that the channels 
Af~*~, A\ , A~^~ and + had contaminations with scattering states of between 20% and 
40%. The A^ + and E ^ channels had bi-torelons contaminations between 20% and 40%. 

Identification of the spin of glueball states on the lattice is non-trivial. Continuum spin 
representations break down into lattice representations at finite lattic spacing. In Tab. 3 
the continuum spin representations are broken down into lattice representations. 



- 8 - 





Figure 6. Glueball masses for the different lattice representations. Open circles refer to results 
obtained in this study, filled squares are quenched results from [3]. 
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Ch. 


^min ^max 


am 


M 2 


X 2 /dof 


mixc 


mixs 


mixT 


A++ 


1-5 


0.839(28) 


0.94(3) 


0.23 


0.6874 


0.2661 


0.0466 


A++* 


1-5 


1.447(55) 


0.98(5) 


0.48 


0.6779 


0.2611 


0.0609 


,1++** 


1-4 


1.75(11) 


0.88(10) 


0.91 


0.6865 


0.2124 


0.1010 


^++*** 


1-5 


1.99(12) 


1.16(14) 


0.91 


0.6727 


0.2055 


0.1217 


At 


1-5 


2.55(39) 


0.99(38) 


0.31 


0.7207 


0.2461 


0.0331 


At 


1-6 


1.30(4) 


0.89(3) 


0.15 


0.3257 


0.0908 


0.5835 


At* 


1-5 


1.57(8) 


0.73(6) 


0.79 


0.4919 


0.1284 


0.3797 


At** 


1-5 


2.10(28) 


0.76(20) 


0.97 


0.6800 


0.1605 


0.1595 


Ar 


1-4 


2.41(47) 


0.81(39) 


0.64 


0.7263 


0.2486 


0.0251 


At 


1-4 


1.94(13) 


0.99(13) 


0.30 


0.7509 


0.1995 


0.0495 


At 


1-4 


2.36(30) 


1.08(33) 


0.06 


0.7086 


0.1708 


0.1206 


At 


1-5 


1.70(10) 


0.76(8) 


0.48 


0.8366 


0.0284 


0.1350 


At* 


1-3 


2.06(10) 


0.79(15) 


0.41 


0.8850 


0.0535 


0.061 


A 2 


1-4 


2.03(20) 


0.66(13) 


0.56 


0.8528 


0.1067 


0.0405 


At* 


1-4 


2.74(56) 


0.97(60) 


1.03 


0.8528 


0.1067 


0.0405 


A 2 - 


1-5 


2.32(25) 


0.98(25) 


1.10 


0.5526 


0.3905 


0.0569 


E++ 


1-5 


1.202(32) 


0.86(3) 


0.97 


0.8117 


0.1622 


0.0261 


E++* 


1-5 


1.606(87) 


0.89(8) 


0.44 


0.8457 


0.1213 


0.0330 


E+- 


1-4 


2.42(40) 


0.98(39) 


0.01 


0.4243 


0.1398 


0.4359 


E-+ 


1-5 


1.42(6) 


0.76(5) 


0.76 


0.4243 


0.1398 


0.4359 




1-3 


1.71(10) 


0.81(8) 


0.57 


0.6943 


0.0143 


0.2914 


E— 


1-5 


1.86(19) 


0.71(13) 


1.1 


0.7984 


0.1586 


0.0430 


E * 


1-5 


2.08(22) 


0.83(18) 


0.92 


0.8411 


0.1464 


0.0126 


T++ 


1-4 


2.03(16) 


0.92(15) 


0.90 


0.9011 


0.0327 


0.0662 


Tt 


1-5 


1.530(71) 


0.86(6) 


0.45 


0.8870 


0.0689 


0.0441 


rpH * 


1-4 


1.65(7) 


0.87(6) 


0.26 


0.9277 


0.0371 


0.0352 


Tt 


1-5 


2.16(26) 


0.89(23) 


0.24 


0.6700 


0.2759 


0.0541 


rri — h* 


1-5 


2.21(28) 


0.85(24) 


1.15 


0.8308 


0.1114 


0.0578 


Tr 


1-5 


2.25(32) 


0.96(31) 


0.33 


0.7504 


0.1645 


0.0851 


Tt 


1-4 


1.238(43) 


0.92(4) 


0.63 


0.9640 


0.0144 


0.0216 


1 2 


1-4 


1.723(70) 


0.92(6) 


0.85 


0.9551 


0.0291 


0.0159 


T +- 
1 2 


1-4 


1.90(16) 


0.83(13) 


0.99 


0.8336 


0.0504 


0.1160 


rp^ * 

1 2 


1-4 


2.28(29) 


0.94(26) 


0.92 


0.9249 


0.0350 


0.0401 


T -+ 
1 2 


1-4 


1.50(8) 


0.76(6) 


0.93 


0.8101 


0.1736 


0.0163 


rri — h* 
1 2 


1-5 


1.69(12) 


0.90(10) 


0.28 


0.3989 


0.5795 


0.0217 


T 2 - 


1-4 


2.35(26) 


1.07(26) 


0.91 


0.5124 


0.3946 


0.0930 



Table 2. Spectrum on the ensemble of Tab. 1. For each state we were able to extract a signal, we 
show the parameters of the fit with Eq. (2.4) and the relative projection of the mass eigenstate on 
the different subsets of operators (mixc rnixg, mixy). Good overlaps |co| 2 ~ 0(1) and x 2 /dof are 
shown. 
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J 




A 2 


E 


Ti 


T 2 





1 














1 











1 





2 








1 





1 


3 





1 





1 


1 


4 


1 





1 


1 


1 



Table 3. Subduced representations J 4- Go of the octahedral group up to J = 4. This table 
illustrates the spin content of the irreducible representations of Qq in terms of the continuum J. 





J PC 


Mass MeV 






Unquenched 


Quenched 






This work 


M&P 


Ky 


Meyer 




o-+ 




2590(40)(130) 


2560(35)(120) 


2250(60)(100) 




2-+ 


3460(320) 


3100(30)(150) 


3040(40) (150) 


2780(50)(130) 




o-+ 


4490(590) 


3640(60)(180) 




3370(150)(150) 




2-+ 








3480(140)(160) 




5-+ 








3942(160) (180) 


o— 


(exotic) 


5166(1000) 










1~ 




3850(50)(190) 


3830(40)(190) 


3240(330)(150) 




2~ 


4590(740) 


3930(40)(190) 


4010(45)(200) 


3660(130)(170) 




2~ 








3.740(200)(170) 




3~ 




4130(90)(200) 


4200(45) (200) 


4330(260) (200) 




1+- 


3270(340) 


2940(30) (140) 


2980(30) (140) 


2670(65)(120) 




3+- 


3850(350) 


3550(40)(170) 


3600(40) (170) 


3270(90)(150) 




3+" 








3630(140)(160) 




(exotic) 




4140(50)(200) 


4230(50) (200) 




o+- 


(exotic) 


5450(830) 


4740(70) (230) 


4780(60) (230) 






5+- 








4110(170)(190) 




0++ 


1795(60) 


1730(50) (80) 


1710(50)(80) 


1475(30)(65) 




2++ 


2620(50) 


2400(25)(120) 


2390(30)(120) 


2150(30)(100) 




0++ 


3760(240) 


2670(180)(130) 




2755(30)(120) 




3++ 




3690(40) (180) 


3670(50)(180) 


3385(90)(150) 




0++ 








3370(100)(150) 




0++ 








3990(210)(180) 




2++ 








2880(100)(130) 




4 ++ 








3640(90)(160) 




6 ++ 








4360(260) (200) 



Table 4. Glueball masses with J PC assignments. The column M&P reports results from Morn- 
ingstar and Peardon [2] from quenched QCD. The column labelled Ky is the data from Chen et 
al. [3]. Meyer's results are from [25]. 
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Meyer and Teper [25, 55] have developed systematic techniques to classify glueball 
masses into spins. Dudek et al. [56] have stressed the importance of correctly identifying 
lattice representations with spin in the charmonium system. Here, we use the simplest 
spin identification. We identity the J = state with the A\, and J = 1 state with the T\ 
representation, and J = 2 states with almost degenerate T2 and E representations. For 
J = 2 and J = 3 states we take a weighted average of the masses in the component lattice 
representations. 

Morningstar and Peardon [2] found glueball states with h (ground and excited), and 
2 *~ quantum numbers. In the PC = — + sector in Fig. 6 there are two almost degenerate 
levels in the E and T2 representations that could come from a J = 2 state. There is one 
level in the A\ representation which is relatively isolated and so is expected to couple to 
h states. There are also degenerate levels in the A2, and T\ representations that could 
be 3 '". This would be interesting, because 3 h is an exotic quantum number, however, 
because there is no signal in T2 representation, this state is not included in the summary 
tables. 

In the PC = sector in Fig. 6, the errors are large because the masses are very 

heavy. Morningstar and Peardon [2] found glueball states with 1 , 2~~, and 3 . Here 
we see evidence for a 2 state with near-degenerate masses in the E and T2 channel. The 
results in Fig. 6 are not consistent with 3 , because there is no signal in the A2 channel. 
There is also a hint of a signal for the spin exotic state in the Ai representation, but 
unfortunately with large errors. 

For the operators with PC = ++ in Fig. 6, there is potentially a rich mixture of 
states. Morningstar and Peardon [2] found glueball states with ++ (ground and excited), 
and 2 ++ . There are two ++ states in the A\ channel. There are two almost degenerate 
states in the E and T2 channel that would correspond to a 2 ++ state in the continuum. 
There is potentially a 4 ++ appearing in the A\, E, Ti, and T2 representations. There is 
also a potential 3 ++ state in the A2 representation, but there is no hint of degenerate states 
in the T\ and T2 representations. 

In quenched QCD, Morningstar and Peardon [2] found glueball states with J PC = 
3 H , 1 H , 2 H , and O" 1 . The summary of the masses in Fig. 6 doesn't show a signal in 
the E~^ , so our results are inconsistent with an exotic 2^ state. There is evidence for 
states with J PC = 3+~, and 0+ - . 

In Tab. 4 we summarise our assignments for J PC quantum numbers to the glueball 
masses from operators, and compare with other results in the literature from quenched 
QCD calculations 2 . In Fig. 7 we plot the summary of the glueball masses. 

To compare with experiment we use the masses in the PDG. However we also include 
the additional masses reported by the Crystal Barrel collaboration, reviewed by Bugg [60]. 
Note that the Crystal Barrel experiment could only find new states in a restricted mass 
range. As explained in [59], some of the results from Anisovich et al. [57-59] are not 
included in the PDG summary tables, because they are neither confirmed nor excluded 

2 We quote results in MeV, but the reader should be aware that different scale settings have been used 
in various cases. For each work, the scale setting is described in detail in the corresponding reference. 
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Figure 7. Summary of the glueball masses compared to experimental meson masses. The 
experimental results are from the PDG [29] and from Anisovich et al. [57-59]. 

by other experiments. We plot the results for 1 = mesons in [57, 58], but plot the 2 ++ 
states from [59]. These additional meson states from Anisovich et al. [57-59] have been 
used to find some experimental glueball candidates [61] from the quenched calculations of 
Morningstar and Peardon. 

Bali [62] has plotted the quenched glueball spectrum with the experimental masses of 
the charmonium system. 

5 Conclusions and future prospects 

The most conservative interpretation of our results is that the masses in terms of lattice 
representations are broadly consistent with results from quenched QCD. We do not see any 
evidence for large unquenching effects, however a definitive calculation requires a continuum 
extrapolation, and the inclusion of fermionic operators. In Tab. 7 we tentatively assign J PC 
quantum numbers to 10 glueballs. 

Of particular note in Tab. 4 is that Meyer and Teper [25] do not see the two spin 
exotic states identified by Morningstar and Peardon [2]. In their summary of the glueball 
spectrum Morningstar and Peardon note that their spin exotic glueball 2" 1 could actually 
be part of 5 H , 7 H , or ll -1 glueball. Mathieu [42] have also compared the results 
for glueball masses from Morningstar and Peardon with those from Teper and Meyer. 
Our result for the mass of the H are consistent with the result from Morningstar and 
Peardon [2], although our errors are large for this heavy state. Meyer and Teper [4, 25] used 
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sophisticated measurement techniques to help assign J quantum numbers to their results 
in terms of lattice representations, but other groups reporting glueball masses have not 
done this. Given the importance of the masses of the oddball glueballs to the experimental 
program of PANDA, this issue needs to be resolved. Future lattice QCD calculations using 
heavy glueball degrees of freedom should use improved techniques to assign quantum 
numbers [4-6]. 

It would be advantageous to lattice QCD calculations if there were no light mesons 
above 3.1 GeV, because this would reduce the complicated mixing of glue and quark degrees 
of freedom. Aesthetically, it would be better to find an isolated glueball, rather than have 
to determine the mixture of glueball and quark degrees of freedom in a state such as the 
_/o(1790). The determination of glueballs with exotic J may be easier for lattice QCD 
calculations, although there is still the problem of the heavy mass and the poor signal to 
noise ratio. 

Unquenched glueball calculations require much higher statistics than for the majority of 
flavour non-singlet lattice QCD calculations. We can use the numbers in Tab. 4 to estimate 
roughly the number of configurations required to achieve a given accuracy. The most 
interesting states in Tab. 4 are the ones with exotic J PC = 0" 1 and 2" 1 quantum numbers. 



Using simple Vy-^configs scann §; we can estimate that we need 550,000 configurations 

to get 50 MeV errors for the spin exotic 0" 1 glueball. 

The MILC collaboration has reported time estimates to produce ensembles of gauge 
configurations with 2+1+1 flavours of sea quarks with the HISQ improved staggered action. 
With the light quark a tenth of the strange quark mass, they estimate the time to generate 
1000 configurations as 35, 128, and 352 Million core hours for the lattice spacings of: 0.09, 
0.06, and 0.045 fm respectively 3 . These numbers suggest that it is still too expensive to 
reduce the errors on the masses of the oddball glueballs. Lattice QCD calculations that 
use anisotropic lattices may be computationally cheaper. 

The heavy glueballs may mix with charmonium states. To study this type of mixing 
will probably require both charm loops in the sea and the computation of disconnected 
diagrams [63, 64] for the charmonium valence states. There are a number of lattice QCD 
calculations, such as the one by the MILC collaboration [65], that include the dynamics of 
the charm quarks in the sea. 

It is an exciting time for lattice QCD, with some unquenched lattice QCD calculations 
done with physical quark masses, and many unquenched calculations resulting in errors at 
the 0(1%) level. The progress on lattice QCD calculations of glueballs is slower. However, 
the start of the PANDA experiment in 2018 provides an important deadline for lattice 
QCD calculations of glueballs. 
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3 Talk presented by Paul Mackenzie at May 2011 review of the lqcd project. 
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